Allometry of the quasi-pipe (qPipe) model for estimating tree leaf area and tree leaf mass applied to plant functional types

The allometry of the pipe model quantifies the approximate proportionality between the tree leaf amount and the stem cross-sectional area at the crown base (ACB). It is useful for estimating and modeling carbon fixation abilities of trees but requires climbing the tree and is thus unsuitable for large-scale studies. Here, we adopted a previously proposed allometry (hereafter the quasi-pipe (qPipe) model allometry) formulating the relationship between the tree leaf amount and a surrogate of ACB, ACB_Est, calculated from tree dimensions measurable from the ground. Using published/unpublished data for 962 trees of 159 species collected between tropical rainforests and boreal forests, we established pipe and qPipe model allometries for evergreen-conifer, deciduous-conifer, evergreen-broadleaf, and deciduous-broadleaf plant functional types (PFTs). For the leaf area per tree (LA), allometric lines on a log–log plane were almost identical among the four PFTs in both models, with slopes of ~ 1. For the leaf mass per tree (LM), however, the allometric lines separated among the four PFTs in both models and had slopes greater than 1, indicating that the proportionality assumed in the pipe model held for LA but not LM. The applicability of the qPipe model in estimating the stand-scale leaf amount was further examined.

www.nature.com/scientificreports/ reason is that the leaf amount of a given tree may not change or even decrease, instead of increasing, despite the DBH and tree height continuing to increase with age 8,9 . This is one reason why both the intercept and slope coefficients of a linear allometric equation change with stand age [10][11][12] .
The allometric equation known as the pipe model of the tree shape 13,14 is suitable for estimating the leaf amount of a tree. The model predicts the tree leaf amount as being approximately directly proportional to the stem cross-sectional area at the crown base (A CB ) or the crown base stem diameter squared (D CB 2 ). As the pipe model uses A CB or D CB 2 as the predictor of the leaf amount, it can express a constant or decreasing leaf amount over time; a tree with a decreasing leaf amount generally has a decreasing crown length because the speed of the rise of crown base location exceeds the height growth owing to the death of lower branches 8,9 . Such a rapid crown-base rise often leads to a decrease in A CB and the leaf amount 8,9 . After the publication of the pipe model, many studies confirmed the validity of the model and related it to the relationship between the tree leaf amount and stem sapwood area at the crown base 15,16 . In this paper, however, we refer to the pipe model as the approximate proportional relationship of the tree leaf amount with the stem cross-sectional area at the crown base A CB or D CB 2 and not the relationship with the sapwood area of a trunk as in other studies. The validity of the pipe model implies that the 'pipe model ratio' 17 (defined as the leaf amount per unit area of the stem cross section or unit sapwood area) is constant irrespective of the tree size. Meanwhile, little attention has been paid to whether the amount of leaves in the pipe model should be expressed on a leaf mass basis or a leaf area basis 16 . The leaf amount in the original pipe model 13 was expressed on a leaf mass basis, whereas some later studies expressed the leaf amount on a leaf area basis. If the pipe model ratio is constant for both the tree leaf area and tree leaf mass, the tree leaf mass per tree leaf area (tree LMA) must be constant between trees of different size in a stand. However, studies of an individual-leaf scale suggest that the LMA of an individual leaf is affected by the light conditions such that it decreases with decreasing light from the canopy surface downward [18][19][20][21] . It is thus expected that the tree LMA is smaller for shorter trees in a stand than for taller trees in the same stand. If so, the pipe model ratio being constant between trees cannot hold for both the leaf area and leaf mass simultaneously. In other words, the proportionality between A CB or D CB 2 and the tree leaf amount cannot hold simultaneously for both the leaf area and leaf mass, unless the tree LMA is constant among trees of different size in a stand. This means that when the allometric line is drawn on a log-log plane, the slope should be different between the allometry for the leaf area and that for the leaf mass.
Care must be taken when estimating the slope of a regression line. In estimating the slope, the use of a regression method called Model II regression is recommended as it calculates reasonable estimates of the slope by considering residuals of both 'x' and 'y' simultaneously 6,22 . Meanwhile, the regression method used for predicting the 'y' value from a given value 'x' is called Model I regression. As the slope of the regression line of Model II regression is known to be greater than that of the Model I regression 22 , allometric analyses for establishing the relationship (Model II) and for prediction (Model I) should be conducted simultaneously.
Another practical problem of the pipe model is that one needs to climb to the base of the tree crown to measure A CB . This issue has partly been solved adopting a method for approximating A CB without the need for climbing 12 . The following relationship was obtained from measurements for 156 trees of Betula ermanii Cham. in three stands 12 : where A BH is the stem cross-sectional area at a 1.3-m height, L C is the crown length (i.e., the difference between the tree height and height of the crown base), and H >1.3 the length of the stem part above a height of 1.3 m. This is equivalent to assuming that there is a virtual paraboloid of revolution that has its apex at the treetop, and cross-sectional areas at the crown base (i.e., at the length L C from the treetop) and at the breast height being A CB and A BH , respectively. It follows from Eq. (1) that A CB can be approximated as 12 where A CB_Est is an estimate of the stem cross-sectional area at the crown base, and R CROWN is the crown ratio defined as L C /H >1.3 . As A BH , L C , and H >1.3 are all measurable from the forest floor, the use of A CB_Est as the surrogate of A CB is practical in estimating the tree leaf amount just like for the pipe model. So far, the application of Eq. (2) has only been reported for two species 12,23 . Both studies showed that the actual measured A CB and its estimate A CB_Est were not statistically different, such that A CB and A CB_Est were treated as though they were convertible predictors of the leaf amount. The allometry for the tree leaf area using A CB_Est as a predictor variable was found to be not statistically different from the pipe model allometry using A CB as a predictor variable 12 . When estimating the LAI of the stand, the tree leaf area for trees whose A CB could not be directly measured in the cited study was estimated by substituting A CB_Est into an existing pipe model allometric equation, namely the A CB -tree leaf area relationship 12 . In contrast, the present study presumes that the pipe model allometry is unavailable owing to a lack of A CB data, and this is why the allometry that uses A CB_Est as a predictor variable is needed. Therefore, we deal with the allometry using A CB_Est and that using A CB (i.e., the pipe model allometry) differently in the present study. Hereafter, the allometry using A CB_Est as a predictor variable for the tree leaf amount 12 is referred to as the quasi-pipe (qPipe) model allometry to indicate the difference. The form of the allometric equations in the present study is given in the caption of Table 1. A comparison between the pipe and qPipe model allometries is presented as equations (S1) and (S5) in Supplementary Information S1 and Supplementary Fig. S1.
The reason why Eq. (2) holds for estimating A CB is not biologically based but is simply based on a general characteristic that the tree trunk tapers toward the treetop (see the Discussion). As Eq. (2) does not require the determination of regression coefficients such as the intercept and slope, if Eq. (2) is found to hold for other  24,25 . Hereafter, the four PFTs are abbreviated as EB, EC, DB, and DC (excluding monocods such as palm-tree species).
In this study, we propose the qPipe model allometry for estimating the stand-scale leaf biomass and LAI in large-scale studies, by pooling the data of the tree species of each PFT. Published/unpublished data for 962 trees collected from tropical rainforests in Southeast Asia, temperate forests of Japan and Tasmania, and boreal forests of Japan and Siberia are used in analyses. The number of trees (species) is 453 (102), 225 (10), 227 (44), and 57 (3) for EB, EC, DB, and DC PFTs, respectively. The list of species and number of trees are given in Supplementary www.nature.com/scientificreports/ Information S3. Hereafter, the term 'leaf mass' refers to the leaf dry mass. The leaf area for broadleaf species is the one-sided area of a leaf. For conifer species, the leaf area in all data used in the analyses is measured as the projected area by detaching all needles from the stem. We first examine (1) whether the relationship between A CB and A CB_EST is 1:1 using standardized major axis (SMA) regression (the smatr library 22 in R 26 ), which is Model II regression suitable for estimating the slope 6 . Then (2) for the four PFTs, the pipe model allometries are explored and compared with qPipe models through SMA regression. Next, (3) to estimate the leaf area and leaf mass and their 95% prediction intervals using Model I regression, the qPipe models are formulated via general linear mixed models (GLMMs), with species being designated as a random effect. This regression method is a Model I regression 6 . In both SMA regression and GLMMs, best models are selected using the Akaike information criteria (AIC). Lastly, (4) to present an example of how LAI estimates for one stand differ between the pipe and qPipe model allometries, a comparison is made between LAI estimates using published 20-year changes (from 21 to 40 years of stand age) in the LAI measured in an evergreen conifer 'hinoki' cypress (Chamaecyparis obtusa (Siebold et Zucc.) Endl.) stand 8 . Using tree census data including data of the crown base stem diameter recorded for each tree each year over the 20 years 9 , the LAIs of each year and their 95% prediction intervals are calculated by adopting three methods; i.e., adopting the pipe model allometry prepared for the stand 8 using the measured A CB as the predictor ('site-specific pipe model'), the qPipe model allometry prepared using A CB_Est as the predictor, where A CB_Est is calculated using the same tree census data ('site-specific qPipe model'; see Supplementary Information S1), and the qPipe model allometry of the PFT obtained in (3) using A CB_Est ('global qPipe model'). In this comparison, the 95% prediction intervals of each LAI estimates are calculated through an error propagation of 95% prediction intervals of the leaf area estimates of individual trees 7,8 . All statistical calculations are performed using R 26 .

Comparison between measured A CB and estimates A CB_Est .
The measured values of A CB (m 2 ) and the estimates A CB_Est (m 2 ) calculated with Eq. (2) were compared among PFTs using SMA regression (Fig. 1). Both A CB and A CB_Est were log-transformed in the regression calculation. The best model was the one having ln(A CB ) and the four PFTs (categorical variable) as the predictor variables ( Table 1). The allometric line was drawn for each PFT. Figure 1 shows that the relationship was almost 1:1 for each PFT. The slopes of the four PFTs ranged between 0.99 and 1.01, without significant differences from 1.0 (Table 1). Pair-wise differences were significant between the elevations (intercepts) of the regression lines of the four PFTs except between the elevations of the EC and DC PFTs and between the elevations of the EB and DB PFTs. However, these differences were small as illustrated in Fig. 1, in which the four regression lines almost overlap. This result suggests that A CB_EST can be used practically as a surrogate of A CB for each PFT.
Comparison between the pipe-and the qPipe-model allometries using Model II (SMA) regression. For the pipe model relationship using SMA regression, the model having ln(A CB ) and PFTs as the predictor variables was selected with the natural-log transformed leaf area (Fig. 2a) and leaf mass (Fig. 2b) as the response variables. For the leaf area (Fig. 2a), the slopes of allometric lines were not significantly different between the four PFTs (likelihood ratio test, P = 0.053), ranging between 1.15 and 1.01 (Table 1). The slopes of broadleaf PFTs (EB, DB) were not significantly different from 1, indicating that the pipe model ratio for the leaf area was constant and   www.nature.com/scientificreports/ the slope follows the pipe model presumption in broadleaf PFTs. In contrast, for conifer PFTs, the difference from 1 was significant in the EC PFT (slope of 1.09, Table 1). The DC PFT had the largest slope (1.15) but the difference from 1 was not significant, probably because only 15 leaf area data were available for the DC PFT, all for Japanese larch (Larix kaempferi (Lamb.) Carr.) trees. Slopes of the conifer PFTs suggest that their pipe model ratio tended to increase with increasing tree size. Pair-wise comparisons among PFTs showed that the difference between PFTs in the elevation of the allometric line was significant only between DC and EC PFTs (Table 1). Despite these differences, regression lines of the four PFTs almost overlapped (Fig. 2a), indicating that pipe model allometries for leaf area were similar across PFTs. In contrast, for the leaf mass, the separation of allometric lines was distinct among the four PFTs (Fig. 2b). Differences were significant in both slopes (likelihood ratio test, P = 0.018) and intercepts (P < 0.000). Elevations of regression lines were greater for conifer (EC, DC) PFTs than for broadleaf (EB, DB) PFTs and were greater for evergreen (EC, EB) PFTs than for deciduous (DC, DB) PFTs [likelihood ratio test for multiple comparison, all P < 0.000; Table 1 (2b)]. Furthermore, the slopes tended to be greater for conifer (EC, DC) PFTs than for broadleaf (EB, DB) PFTs (Table 1 (2b)). More importantly, all four PFTs had significant slopes greater than 1, indicating that the pipe model ratio for the leaf mass differs by tree size, with larger trees having a greater pipe model ratio.
The qPipe models that used ln(A CB_EST ) and PFTs as the predictors (Fig. 2c, d) gave results similar to the pipe model (i.e., Fig. 2a, b) with a few exceptions: for leaf area (Fig. 2c), the difference from 1 in the slope was not   Table 1 for coefficients of the power function of each allometry. The significance of the difference in elevations of the regression lines between the four PFTs are given in Table 1 2c)). Although differences in elevation were significant between EC and DB PFTs (likelihood ratio test for multiple comparison, P < 0.000) and between EC and DC PFTs (P = 0.027), the four allometric lines overlapped around the center of the data distribution (Fig. 2c) as in the case of the A CB -leaf area relationship (Fig. 2a).
For the leaf mass, the qPipe model had clear differences in regression lines (Fig. 2d), as in the case of the pipe model (Fig. 2b). Elevations were significantly different among PFTs (likelihood ratio test, P < 0.000) with conifer PFTs (EC, DC) having higher elevations than broadleaf (EB, DB) PFTs (likelihood ratio test for multiple comparison, P < 0.000) as in the pipe model (Table 1(2d)). Additionally, the slopes tended to be greater for conifer (EC, DC) PFTs than for broadleaf (EB, DB) PFTs and were significantly different from 1 except for the DB PFT (Table 1(2d)).
In summary, for the leaf area, both the pipe model and qPipe gave regression slopes closer to 1, with little apparent difference in the allometric line among PFTs. For the leaf mass, in contrast, slopes of the four PFTs were greater than 1 in both the pipe model and the qPipe model, and their allometric lines tended to differ in elevation.
Prediction of the tree leaf area using the qPipe model allometry. To predict the tree leaf area, the allometric relationship of the qPipe model was obtained through model selection with GLMMs (Fig. 3b). For comparison, the results of the pipe model using GLMMs (Fig. 3a) are also presented. In both cases, the PFT was not selected as a predictor of the Model I regression, in contrast with the results of Model II regression (Fig. 2). Selected models were those with ln(A CB ) or ln(A CB_EST ) as the fixed effect variable and species as the random intercept term (Fig. 3a, b, Table 1). This indicates that, for predicting leaf area using A CB or A CB_EST , differences in the PFT are negligible as far as our collective data sets are concerned. In other words, the qPipe model (Fig. 3b) is applicable to the prediction of the leaf area of a tree without the need to consider differences in the PFT. Note, however, that different species would have different elevations of the allometric line as the species was selected as the random effect on the elevation. Ninety-five percent prediction intervals are drawn in Fig. 3 using an approximating equation, which is given in Table 2 (see also Supplementary Information S2). The intervals are seen to be slightly wider for the qPipe model (Fig. 3b) than for the pipe model (Fig. 3a).
Prediction of the tree leaf mass using the qPipe models. The results presented in Fig. 2b, d clearly show that, in predicting the tree leaf mass, the four PFTs have different qPipe model allometry equations. Hence, the regression calculation was performed for each PFT independently, allowing the calculation of prediction intervals. Models with species as the random intercept term were selected for the EC, EB, and DB PFTs (Fig. 4a, c, d) whereas ordinary least-squares (LS) linear regression without a random effect was selected for the DC PFT (Fig. 4b, Table 1). Consistent with the results of Model II regression (Fig. 2d), the elevation of the allometric line tended to be greater for the evergreen trees than for the deciduous trees in both conifer and broadleaf PFTs (Table 1).    Table 1 for coefficients of each allometric equation, Table 2 for the coefficients of the equation approximating the 95% prediction intervals, and Supplementary Information S1 and S2 for the method of approximating the prediction intervals equations.  Fig. 4b, the upper and lower 95% prediction intervals of an estimate of LA* or LM* (= Y ) for a given ln(A CB *) or ln(A CB_Est *) (= X) are given by exp[ln(Y) ± ΔPI 95% ], where ΔPI 95% = (C 0 + C 1 × X + C 2 × X 2 ). For Fig. 4b, ΔPI 95% for ln(Y) is calculated with the formula of the 95% prediction interval in ordinary leastsquares regression; i.e., ΔPI 95% = {C 0 × C 1 × sqrt[1 + 1/C 2 + (ln(X) − C 3 ) 2 /C 4 ]}.    Table 1 for coefficients of each allometric equation, Table 2 for the coefficients of the equation approximating the 95% prediction intervals, and Supplementary Information S1 and S2 for the method of determining/approximating the prediction interval equations. www.nature.com/scientificreports/ Example of the comparison of the pipe and qPipe models for the LAI. Lastly, the three methods of stand LAI estimation were compared. A 20-year change in the LAI of an evergreen conifer hinoki (C. obtusa) stand was estimated using the measured A CB of each tree and the original 8 site-specific pipe model allometry prepared for the particular stand, the calculated A CB_Est and the site-specific qPipe model allometry, and the calculated A CB_Est of each tree and the global qPipe model allometry for the leaf area (shown in Fig. 3b). The LAI estimated with the site-specific qPipe model showed interannual changes similar to those of the LAI estimated using the site-specific pipe model, though there were differences (but not significant) in the LAI (Fig. 5a). A positive relationship between the LAI and 6-year running average of the mean summer temperature, which has been reported previously 8 , was reproduced with the site-specific qPipe allometry (Fig. 5b). However, the LAI estimated with the global qPipe allometry could not reproduce the relationship (Fig. 5a, b), especially at earlier stand ages. These results demonstrate that qPipe allometry is useful if the site-specific qPipe allometry is available together with measurements of the tree height, crown-base height, and DBH for each tree, whereas the global qPipe allometry may need to be carefully considered when applied to the LAI estimation of an unknown stand, as we discuss later.

Discussion
Leaf area versus leaf mass. Although the main objective of the present study is to examine the applicability of the qPipe model, the study also advances our understanding of the pipe model. An important presumption of the original pipe model is that the amount of leaves per unit cross-sectional area of the trunk at the crown base, or the pipe model ratio, is constant 13 . We show that this appears to hold for the leaf area but not for the leaf mass, as long as the stem cross-sectional area A CB is considered as a predictor variable. In the Model II regression analysis (Fig. 2a, Table 1), the slope of approximately 1 for PFTs except for the EC PFT indicates that the pipe model ratio for the leaf area (leaf area per unit A CB ) was statistically constant irrespective of A CB . In contrast, for the leaf mass, the slope of the Model II regression was significantly greater than 1 for each PFT (Fig. 2b, Table 1). This indicates that the pipe model ratio for the leaf mass (leaf mass per unit A CB ) increased with increasing A CB . The difference in the slope of a PFT between the leaf area and leaf mass can be ascribed to a tendency that taller trees tend to have a greater tree LMA than shorter trees, as mentioned in the Introduction. Figure 2a also shows an interesting result for the leaf area that there was little difference in the pipe model allometries, if at all, between PFTs. This implies that the pipe model ratio for the leaf area remained approximately constant irrespective of the PFT as long as A CB is used as the denominator of the pipe model ratio. However, the use of the sapwood area in the stem cross section, rather than the stem cross-sectional area A CB , has been regarded as a reason why the pipe model holds, as the former can biologically explain the hydraulic functioning 15,16 . Additionally, the xylem structure for the water transport system is known to vary between the four PFTs 27,28 . For example, the gymnosperm (EC and DC) species depend on tracheids for water transport and do not have Each of the three allometries takes the form LA = CF × Elevation × X Slope , where LA (m 2 ) is the leaf area per tree, X represents A CB (the site-specific pipe model) or A CB_Est (the site-specific and global qPipe models) of a tree, and CF is the factor correcting for bias regarding log-transformed regression. The values of the CF, elevation and slope differ between the three methods. See Supplementary Fig. S1 and Supplementary Table S1 in Supplementary Information S1 for coefficients of these allometries. For the site-specific pipe and qPipe models, allometric equations were formulated using 47 trees sampled near the stand where the LAI was estimated 8 . The allometry of the global qPipe model is the same as in Fig. 3b and Table 1(3b). LA for each tree each year in the hinoki stand was calculated by substituting measurements of A CB or A CB_Est of each tree each year into the three different allometries. The LAI for each year was then calculated by taking their sum and dividing by the stand area. (b) Relationship between the LAI of each year and the moving average of the mean temperature in July and August for the last 6 years (T JA6y ), where the LAI was calculated using the site-specific qPipe allometry (open squres) or the global qPipe allometry of Fig. 3b (open circles). In the original paper, a significantly positive relationship was obtained 8 . For coefficients of the regression, see Supplementary Table S1 in Supplementary Information S1. www.nature.com/scientificreports/ vessels, unlike angiosperm (EA and DA) species. Even between angiosperm species, the vessel structure of the xylem on the stem cross-section varies; e.g., between ring-porous wood (e.g., some deciduous oak spp.), diffuseporous wood (some beech spp.), and radial-porous wood (some evergreen oak spp.). Considering these possible variations in xylem hydraulic properties, it seems curious that the slope of approximately 1 held for the leaf area-A CB relationship with almost no difference among the four PFTs. We do not have an explanation for this as yet. Considering that the leaf area rather than the leaf mass is likely the primary quantity of leaves relating to photosynthesis 20,29 , trees may adjust to maximize water conduction efficiency for photosynthesis and transpiration irrespective of the xylem structure 30 , which may have resulted in a convergence to a similar proportionality among PFTs in the leaf area-A CB relationship. For the pipe model relationship to hold, stem thickening growth must occur in such a way that the stem cross-sectional area or sapwood area is coordinated to be proportional to the leaf amount. However, it is not known what types of signal are involved in this coordination. One possibility would be the effect of the plant hormone auxin (indole-3-acetic acid, IAA), which plays an important role in determining the radial growth of a stem by affecting the activity of cell division in the cambial meristem 31 . As IAA is produced in apical buds and growing shoots in the crown and transported to the cambial region of the stem, the IAA level in the cambial region is known to be affected by the crown size [32][33][34][35] . This may be a cue to the connection between the leaf amount and the stem cross-sectional area. Meanwhile, a twig-scale study suggested a mechanism that allows the pipe model relationship to hold in terms of the source-sink relationship 36 . Further study combining several biological mechanisms may be necessary to solve the abovementioned questions.
How and when does the qPipe model work? Figures 1 and 2 show that the qPipe model provided allometry comparable to that obtained using the pipe model. However, whether the qPipe model can fairly estimate the actual A CB or not depends on the actual form of the tree stem. Here, we explain this using the stem taper equation, a mathematical expression of the change in stem thickness from breast height to treetop. In the field of forest sciences, the shape of a tree stem is often expressed by a simple equation 37 . The relative taper of the stem between the crown base and breast height is expressed as where L is the distance from the treetop along the stem (relative scale; 0 < L ≤ 1, where L = 1 at breast height), A is the stem cross-sectional area at L (relative scale; 0 < A ≤ 1, where A = 1 at breast height), and X is a scaling parameter indicating the degree of stem taper. The stem shapes of a cone, paraboloid of revolution, and cubic solid of revolution correspond to X = 1, 2, and 3, respectively (curved lines in Fig. 6a). As mentioned in the Introduction, the calculation of A CB_Est in the qPipe model is based on the assumption of a paraboloid of revolution; i.e., X = 2 in Eq. (3). It follows that, if the qPipe model is applied to a tree stem having X < 2, then A CB_Est > A CB (i.e., an overestimation, as L 2/2 > L 2/X ), and if applied to trees having X > 2, then A CB_Est < A CB (i.e., an underestimation,  8,9 were used in the analysis. (a) Relative crown length and A CB relative to A BH for each live tree at 21 (filled circles) and 40 (open squares) years of stand age in an even-aged population of a hinoki cypress stand 8,9 . X = 1, 2 and 3 corresponds to the stem shape of a cone, paraboloid of revolution (Eq. (2)), and cubic solid of revolution, respectively. As the line of X = 2 corresponds to the qPipe model, data points on the line of X = 2 imply that A CB and A CB_Est estimated by the qPipe model are in agreement. Pink crosses represent data of trees that died between 21 and 40 years of age; the marks are the data 1 year before their death. Such trees are found to be overestimated. (b) Relative error of the stem cross-sectional area at the crown base (= A CB_Est /A CB ). www.nature.com/scientificreports/ as L 2/2 < L 2/X ) (Fig. 6a). Hence, concordance between A CB and A CB_Est (Fig. 1) depends on the stem shape that a tree actually had, or the X value in Eq. (3). This is a possible source of estimation error of A CB_Est for an individual tree. However, in estimating the total leaf amount of a stand, what should be considered further is the balance between overestimation and underestimation among the trees in the stand. The stand-level leaf amount may be overestimated if there are many trees with X < 2 and underestimated if there are many trees with X > 2. The data points of Fig. 6a illustrate an example of actually measured A CB for all trees in an even-aged evergreen conifer stand at 21 and 40 years of stand age 8,9 . At 21 years of stand age, a majority of the trees had A (≡ A CB /A BH ) at the crown base in the region X < 2 (Fig. 6a), indicating that their A CB_Est was an overestimation of the actual A CB (Fig. 6b). This tendency changed at 40 years of stand age, when there were still some overestimated trees but many underestimated trees also (Fig. 6a, b). This change is ascribed to a change in the stem taper with tree growth under competition. The good estimation of the LAI by the site-specific qPipe model with non-significant differences in the LAI (Fig. 5a) was partly due to the overestimation and underestimation offsetting each other between trees. It is likely that trees of a stand vary in terms of X in Eq. (3).
A study of four conifer species in North America 37 reported that there were more trees with X = 2 than those with X = 3. This may suggest that the assumption of X = 2 in the qPipe model is likely a moderate assumption for estimating the stand LAI due to the presence of trees with both X < 2 and X > 2 in a stand. Additionally, the qPipe model provides over/underestimations due to other factors. For example, many long-lived 'giant' trees 38,39 are known to undergo treetop breakage because of natural phenomena such as lightning strikes and strong winds. Such trees are expected to have a less tapered, cylindrical stem (i.e., X > 2 in Eq. (3)) owing to the loss of their top part and would have a much higher H value if breakage had not occurred. The leaf amount for such top-broken trees may result in underestimation with the qPipe model. Another situation that may lead to error is a species having leaf turnover without turning over branches; evergreen conifers such as black spruce (Picea mariana (Mill.) B.S.P.) turnover needled shoots with epicormic branching (re-sprouting of new shoots) on a primary branch whereas older shoots on the same branch wither and fall 40,41 . For such a tree, the amount of leaves per tree may be in a steady state, whereas A CB may increase annually without a change in its crown-base location. If this happens, the leaf amount per unit A CB should decrease with age. These problems apply to not only the qPipe model but also the pipe model. The applicability of the qPipe model and pipe model needs to be carefully considered for trees having an atypical crown structure.

Conclusion
With our collective dataset, the constant pipe model ratio assumed in the pipe model held for LA but not for LM, which likely reflects changes in the tree LMA with tree size. We proposed qPipe model allometry considering the difference in PFTs, where measurements of the tree height, crown-base height, and DBH are used in estimating the leaf area and leaf mass of an individual tree. As these measurements can be made in the field and the 95% confidence intervals of the estimated individual leaf area or mass are calculable, the stand LAI and its confidence intervals can be estimated for a stand for which the stand-specific allometry of the leaf area or leaf mass is not available, as is often the case in large-scale studies. Although data for ascertaining the applicability of the qPipe model to LAI estimation are still limited, we believe that the investigation of the applicability of the present method is worth continuing in future work.

Materials and methods
Data sources. Many of the data of the present study are taken from the BAAD dataset 42 , which is a metadataset comprising many published/unpublished data of individual trees. However, while choosing data suitable for the present study, we found that a part of the BAAD data had been incorrectly copied from original data. We referred to the original data source in adding/correcting such data. Additionally, other data sources were collected and combined with the BAAD data.
With the combined data, we selected individual trees that satisfied the conditions of a tree height greater than 1.3 m, breast height of at least 1.3 m, crown-base height of at least breast height, either the stem diameter or cross-sectional area being recorded for both the breast height and crown-base height, and either the leaf area (LA, m 2 ) or leaf dry mass (LM, kg) being recorded for an individual tree. In the process of regression calculation, we carried out Grubb's test (R package "outliers" 43 ) to remove outliers from each of the ln(A CB )-ln(LA) and ln(A CB )-ln(LM) relationships until no outliers remained. If necessary, Grubbs' test was conducted for each PFT or each species. Data for palm trees (monocots) were excluded. The list of all of the data used in the analyses of the present study is given as the Supplementary file "AllData.csv", the explanation of which is given in Supplementary Information S4. Data analyses. All statistical calculations were made using R 26 . For Model II regression (Figs. 1, 2), the SMA regression function available from the package smatr 44 in R was used, where the "robust" setting was adopted to lessen the effects of outliers. Model selection was performed using the AIC. As the current version of smatr does not allow the calculation of random effects, PFTs were treated as the predictor by pooling all species in each PFT. In Model I regression used for predictions (Figs. 3, 4), ordinary least-squares or the mixed models package lme4 45 was used, both of which allow the calculation of prediction intervals. In GLMMs with log-transformed variables, a Gaussian error distribution was chosen for normality. Models with the PFT as either a fixed effect or random effect were compared to select the best models using the AIC. The species in each PFT were treated as a random effect, where model selection was made by comparing models including species only in a random-effect intercept term with those including species in both intercept and slope terms. Some models with error messages (e.g., the model failed to converge) were discarded irrespective of the AIC. As a result, in all mixed models of the present study (Figs. 4, 5a, c, d), the species were treated as the only intercept term of the random effect. www.nature.com/scientificreports/ Note that the prediction intervals for the mixed models in Figs. 3 and 4 were first calculated with lme4. The prediction intervals were then approximated with the quadratic equation function of A CB or A CB_Est (Table 2) to allow the calculation of prediction intervals for other studies. This approximation was performed using the nonlinear regression method nls() in R. An example of the results of calculation is given in Supplementary Fig. S2.
As the regression calculation for Model I regression (Figs. 3, 4) was made through natural-log transforming the predictor and the responses of a power function, the bias arising from the log-transformation was corrected by multiplying the power function with a correction factor CF 12,46 . The CF values are given in Table 1. See also equation (S5) in Supplementary Information S1 for an example of the equation.